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Magnetic susceptibility, magnetization, specific heat and positive muon spin relaxation (/iSR) 
measurements have been used to characterize the magnetic ground-state of the spinel compound 
CuGa204. We observe a spin-glass transition of the S=f/2 Cu'^^ spins below Tf = 2.5K charac- 
terized by a cusp in the susceptibility curve which is suppressed when a magnetic field is applied. 
We show that the magnetization of CuGa2 04 depends on the magnetic history of the sample. Well 
below Tf , the muon signal resembles the dynamical Kubo-Toyabe expression refiecting that the spin 
freezing process in CuGa2 04 results in a Gaussian distribution of the magnetic moments. By means 
of Monte-Carlo simulations, we obtain the relevant exchange integrals between the Cu^^ spins in 
this compound. 
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I. INTRODUCTION 

Although spin glasses have been extensively studied 
in the past years, there is still no consensus about the 
ground-state and dynamics in these systems (for an in- 
troduction see e.g. K.H. Fisher and J.A. Hertz [ |). It 
is generally accepted that both site-disorder and com- 
petition between the magnetic moments are necessary 
to produce a low-temperature siate where the spins are 
frozen along arbitrary directional. Examples of such sys- 
tems are metallic spin-glasses where magnetic impurities 
are randomly diluted in a noble metalo. For this par- 
ticular class of materials competition between the mag- 
netic moments is the result of the jIludermann-Kittel- 
Kasuya-Yosida (RKKY) interactiorD where ferromag- 
netic and antiferromagnetic exchange interactions alter- 
nate as a function of distance between neighboring spins. 
The RKKY interaction cannot be invoked for localized 
magnets and the spin-glass transition in these systems 
must be realized by other mechanisms. Typical insulat- 
ing spin-glasses of this kind are the alloys Eua;Sri_a;S. 
In the a;=0-limit, EuS is a well-known example of an 
isotropic 3-dimensional Heisenberg ferromagnet. The ex- 
change integrals have been determined by inelastic neu- 
tron scattering in this material with the result that fer- 
romagnetic nearest-neighbor exchange interaction|-|Com- 
petes with next-nearest antiferromagnetic couplingB. Di- 
luting non-magnetic Sr for Eu ensures bond-randomness 
and the conditions for obtaining a spin-glass state- are 
fulfilled in a large range of impurity concentrationsQ, in 



qualitative agreement with the molecular-field theory of 
Edwards and Andersonfl. De Seze pointed out that a 
spin-glass phase transition can occur in a geometrically 
frustrated systenx.with Ising spins and antiferromagnetic; 
interactions onlyO. Following De Seze's work, VillainEI 
proposed that spin glasses can be obtained in materi- 
als with geometric frustration and Heisenberg-type ex- 
change interactions like cubic spinels. These compounds 
have the chemical formula AB2O4. The chemical struc- 
ture of spinels consists of both tetrahedral and octahe- 
dral sites. The number of crystallographic sites is larger 
than the number of A and B cations in the chemical for- 
mula, so that the cations generally distribute randomly 
among the available atomic positions. In particular, this 
random distribution of cations determines in a large ex- 
tent the nuerowave relaxation properties of the spinel 
compoundal2l. When both sublattices are occupied by 
magnetic ions the ground state is a ferrimagnet. The B 
sublattice builds connected tetrahedra and antiferromag- 
netic interactions induce topological frustrationstll which 
can lead to a spin-slass state when non-magnetic impuri- 
ties are introducedEI. The dominant magnetic interaction 
in most of these materials is antiferromagnetic and con- 
nects spins between the A and B sublattices while the 
A-A and B-B exchange interactions are comparatively 
small. However, intersublattice exchange constants can 
modify the magnetic phase diagram originally calculated 
by Villain and for real systems the situation is usually 
complicatedllj. Althou^j-, spin-glass transition has been 
found in diluted spinelaij, spin-glass transition in pure 
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cubic spinels is less common. 

In this article, we report magnetic susceptibility, mag- 
netization measurements in fields up to 50 KOe, spe- 
cific heat and muon-spin relaxation (/iSR) measurements 
in the cubic spinel CuGa204. The results show that 
CuGa204 undergoes a paramagnetic to spin- glass phase- 
transition at Tf = 2.5K. By means of Monte-Carlo sim- 
ulations, the relevant exchange interactions are obtained 
for CuGa204. We show that the formation of a spin- 
glass ground-state in CuGa204 is due to the Jahn- Teller 
character of the Cu^+ ions. Specifically, in a field of 
octahedral symmetry, the Jahn- Teller effect distorts the 
electronic d levels of the Cu^+ which become split by the 
effect of the crystal field into a threefold degenerate level 
and a twofold degenerate one. In such compounds there 
is an interaction between the electronic system with the 
underlying lattice which very often leads to a structural 
phase transition. Typical compounds exhibiting cooper- 
ative Jahn- Teller distortion are found in e.g. perovskites 
(KCuFg, LaMnOa), spinels (CuFe204, Mn304), rutiles 
(CrF2, CUF2) or garnets (Ca3Fe2Ge30i2). The struc- 
tural phase transition can be accompanied by orbital- 
ordering of the d-electrons which in turn influences the 
nature of the exchange interaction. The important role 
of the Jahn- Teller effect in forming the magnetic ground- 
state in the perovskite manganites which exhibit colossal 
magnetoresistance (e.g. see_[14 ] and references therein) 
and in cuprates {e.g. [ |l^,|l^ ] and references therein) is 
currently a subject of intense investigation both theoret- 
ically and experimentally. In that respect, we note that 
the influence of the Jahn- Teller effect on the properties 
of the magnetic, insulators is discussed in detail by Kugel 
and Khomskyll3. 



II. EXPERIMENTAL DETAILS 

A. Sample Preparation 

Single crystals of CuGa204 were grown by spontaneous 
crystallization starting from a CuO— Ga203 solution melt 
in PbO - O.64B2O3 - 0.5Na2O. After slowly coohng the 
melt to room temperature, single crystals of typical size 
3x3x3 mm'^ and of octahedral shape were obtained. X- 
Ray diffraction analysis showed that the CuGa204 crys- 
tals used for the present experiments are cubic spinels 
with both copper and gallium ions randomly distributed 
in the A and B sublattices in agreement with previous 
diffraction investigation (see Ref. [ |l^). The chemical 
structure of CuGa204 is described by the space- group 
Fdim with lattice constants a=8.39A at room tempera- 
ture. 



B. Magnetic Measurements 

The magnetic susceptibility and magnetization mea- 
surements were performed with a commercial MPMS 
Quantum Design SQUID magnetometer together with 
an AC-susceptibility option at ICMA, Spain. The am- 
plitude of the AG-magnetic field was set to 4.5 Oe with 
the frequency of the field being varied between 0.1 Hz and 
990 Hz. The measurements were carried out in the tem- 
perature range 1.7-300 K and in applied magnetic fields 
up to 50 kOe. Additional measurements of the magnetic 
susceptibility in the temperature range T=4.2K to 120K 
were performed at the Institute of Physics, Krasnoyarsk, 
using a home-built SQUID magnetometer. 



C. Specific Heat Measurements 

The specific heat measurements were performed with a 
commercial PPMS device (Quantum Design) in the tem- 
perature range 1.8 K< T < 10 K. We used a small single 
crystal of mass ~4.15mg. The raw data were corrected 
for the copper host and glue, which were measured sepa- 
rately. We did not attempt to subtract the phonon con- 
tribution, as it is expected to be small at low tempera- 
tures. 



D. Muon-Spin Relaxation 

The /zSR experiments were performed on the LTF 
spectrometer at the Paul-Scherrer Institute, Switzerland. 
The data were recorded using the zero-field method 
which is very sensitive to determine both static and dy- 
namic effects in spin-glasseso. Additional measurements 
were performed as a function of applied magnetic field. 
In that case, the sample was zero-field cooled. The sam- 
ple we used for the present experiment consists of about 
50 pieces of the above described crystals which were 
glued on a silver plate. The sample was enclosed in a 
top-loading "^He-^He dilution cryostat and the measure- 
ments were carried out in the temperature range 650 mK 
< T < 10 K. 



III. MAGNETIZATION, SUSCEPTIBILITY AND 
SPECIFIC HEAT RESULTS 

Figure|l| shows the result of the magnetic susceptibil- 
ity measurements with an AC-frequency of 19 Hz and 
an excitation amplitude of H=4.5 Oe. For temperatures 
higher than T=20 K, the magnetic susceptibility is well 
reproduced by the Gurie- Weiss law. Upon lowering the 
temperature below T=20 K, the magnetic susceptibility 
increases continuously. The real part of the magnetic 
susceptibility shows a cusp at Tf ~ 2.5K which is inde- 
pendent of the relative orientation of the magnetic field 
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with respect to the crystal axes. The imaginary part of 
the magnetic susceptibility also exhibits a maximum at 
the same temperature. To understand the nature of the 
maxima appearing in the susceptibility curves, the mag- 
netization in CuGa204 was determined as a function of 
applied magnetic field and for different magnetic histories 
of the samples. As an example, Fig. ^ shows the magne- 
tization curves obtained in CuGa204 after zero-field and 
field-cooling, respectively. For the latter case, the sample 
was cooled in a magnetic field of H=100 Oe applied along 
the [001] crystal axis. It is evident from the figure that for 
temperatures below Tf ~ 2.5K the magnetization shows 
a temperature hysteresis, which is an usual characteristic 
for the formation of a spin-glass state. The results of the 
magnetic susceptibility measurements taken for different 
magnetic fields are presented in Fig. ^ which shows that 
magnetic fields larger than H=:5 kOe suppress the cusp 
observed at Tf in zero-magnetic field. Figure ^ shows 
that the increase of the magnetic moment as a function 
of magnetic field is far from saturation at the maximum 
field of 5T. In Fig. || the temperature dependence of 
the transition temperature Tf of the spin-glass transi- 
tion is observed to increase as a function of increasing 
AC-frequency. The above experimental results all indi- 
cate that the Cu^+ magnetic moments in CuGa204 un- 
dergo a phase transition to a spin-glass ground state be- 
low Tf ~ 2.5K. This is also confirmed by the calorimetric 
measurements performed in zero-magnetic field for this 
compound. A plot of the specific heat Cp/T in CuGa204 
is shown in Fig. ^. The data do not show any indica- 
tion of a phase transition to a 3-dimensional ferro- or 
antiferro-magnetic ordered state. However, a broad max- 
imum is observed around T = 2.5 K followed by a slow 
decay toward high temperatures. This particular behav- 
ior of the specific heat as a functioH-|Of temperature is 
reminiscent of a spin-glass transitionta. 



IV. DISCUSSION OF THE BULK 
MEASUREMENTS 

A spin-glass state is characterized by an assembly of 
magnetic moments which are frozen along random and 
arbitrary directions in space below a specific transition 
temperature Tf. Because of the non-ergodicity of the 
system, the phenomenon is irreversible. The macroscopic 
magnetization of a spin-glass system is equal to zero in 
the absence of a magnetic field. On the other hand, cool- 
ing a spin-glass in an external magnetic field transfers the 
system into a metastable state with a non-zero magne- 
tization value. For temperatures above the spin-freezing 
temperature Tf , the magnetic moments are in a paramag- 
netic state and consequently the temperature dependence 
of the magnetic susceptibility follows the Curie- Weiss law 



where C = Ng'^ii^S{S + l)/3fcs is the Curie constant 
and 9 the paramagnetic Curie temperature. N is the 
magnetic moment density, g the Lande factor and /is 
the Bohr magneton. S corresponds to the spin value of 
Cu^+ and fc^ is the Boltzmann constant. From the mag- 
netic susceptibility measurements presented above, we 
obtain for CuGa204 the values 1/C=0.34 emu K/mol, 

9 = —8 and fXcS = g-\/S(S + l)/iB = 1.65/iB. The mag- 
netic susceptibility is related to the Edwards- Anderson 
(EA) parameteiuq = limt,^oo[(Si(t)Si(0))]av. through the 
relation 
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C 
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T-9{l-q{T)y 

According to the percolation theory of Kirkpatricl0, the 
EA-parameter follows a power law q(T) cx (1 — T/Tf)'' 
close to the spin-glass temperature Tf with [3 equal to 
0.39. However, near Tf, we found the value /? = 0.16 
in CuGa204. The frequency dependence of the spin 
freezing temperature Tf is a characteristic feature of 
the spin-glass state. It has been experimentally found 
that in spin-glasses Tf decreases with increasing AC- 
frequency. A quantitative measure of the frequency shift 
is obtained from (ATf/Tf)/Alog(cj)=0.067. This is very 
similar to other insulating spin-glasses such as (EuSr)S 
or (FeMg)Cl2 [ ^ and one order of magnitude lower 
than in canonical metallic spin glasses, or one order of 
magnitude larger than in a superparamagnet. The mag- 
netic field and temperature dependence of the magneti- 
zation M(H,T) for a spin-glass system with Heisenberg 
spins of dimension m has been calculated by Toulouse 
and Gabaucd. Within the Sherrigton-Kirkpatrick (SK) 
model and for temperatures below Tf , the expression for 
the magnetization is accordingly given by 

M/ff = 1- [4/(to + 2)]1/33/i'*/V4 (3) 

where h = g^sH /ksTf is the reduced magnetic field. 
Figure ^ shows a comparison between the experimental 
measurements obtained in CuGa204 and Eq.|^. There is 
a good agreement between theory and experiment for re- 
duced magnetic fields below h = 0.6. For larger values of 
the reduced magnetic field h, however, there is a signif- 
icant discrepancy between the observed and calculated 
values of the magnetization, which can be attributed to 
the fact that the SK theory is based on the mean-field 
approximation and makes use of infinite long-range inter- 
actions between the spins. On the contrary, the magnetic 
exchange interactions in CuGa204 have short-range char- 
acter. Moreover, the domain of validity of the mean-field 
theory for a spin-1/2 compound is unclear, in particular 
close to the transition temperature Tf where the critical 
fiuctuations become important. In the same spirit, the 
freezing temperature Tf shows a much more pronounced 
magnetic-field dependence than predicted by mean-field 
theory 

/,2/3 ^ [4/(^ + 2)]i/3(l - Tf{H)/Tfm, (4) 
as evidenced in Fig.Bl 
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V. fiSR RESULTS 

To get more insight into both the static and dynamic 
properties of the Cu^+ magnetic moments in CuGa204, 
we have measured the muon-spin relaxation above and 
below the freezing temperature Tf in this material. Gen- 
erally, the magnetic interactions probed by the implanted 
spin-polarized muon are detected by monitoring the 
asymmetric emission of positrons arising from the weak 
decay of the muon. Recording the positron rate N{t) as 
a function of muon life-time yields 



N{t) = N{0) exp(-t/r)[l + AG^t)] 



(5) 



where A is the initial muon asymmetry parameter. The 
product AGz{t) is often called the /iSR signal. In addi- 
tion, the function Gz (t) can be associated with the muon 
spin auto-corrclation fmiction, i.e. 



(S(t)S(O)) 

52 (0) ■ 



(6) 



where S is the spin of the muon. Typical zero-field /iSR 
signals measured in CuGa204 are shown in Fig. ^. Above 
T ~ 3.8 K (and at least up to 10 K), the data are best 
described by assuming for Gz (t) the form 



Gz,para{t) — Gxrit) ' Ges{t) 



(7) 



with G/fTr-cepresenting the familiar Kubo-Toyabe (KT) 
expression^ 

GKT{t) - ^ + ^(1 - ALi')exp(-iALi') , (8) 

and with Ges given by 

Ge.(i) = exp[-(At)''] . (9) 

The form of Gz.para{t) points for the occurrence of two in- 
dependent channels of depolarization acting on the muon 
spin. The first channel, giving rise to the KT function 
Gkt, originates from the nuclear dipolc moments (Ga 
and Cu isotopes). The internal fields of this contribution 
are assumed to be Gaussian distributed in their values, 
randomly oriented and static within the /j,SR time win- 
dow. The parameter A'^^/j'^ represents the second mo- 
ment of this field distribution due to the nuclear spins 
along one Cartesian axis (7^ = 27rT3. 553879 kHz/G is 
the gyro-magnetic ratio of the muon). The second chan- 
nel, described by the function Ges(i), which will be dis- 
cussed in details below, represents the contribution aris- 
ing from the fluctuating electronic Cu spins. 

At T = lOK the muon dcpolarisation can be satisfacto- 
rily described by assuming Ges{t) = 1, i-e. Gz.para{t) = 
GKT{t) with Ans = 0.16(1) MHz (see also Fig. |), indi- 
cating that the fluctuations of the electronic spins are still 
too fast to be observed in the ^SR time window. How- 
ever, upon cooling the sample below T = 10 K and down 



to 3.8 K, the fluctuation rate of the electronic spins de- 
creases and the muon-spin depolarisation becomes grad- 
ually dominated by the Gesit) contribution. Figure |l^ 
represents the temperature evolution of the depolarisa- 
tion rate A. Whereas the exponent (3 remains constant 
in this temperature interval (i.e. (3 ~ 0.78), the depo- 
larisation rate exhibits a marked critical-like divergence, 
which must be taken as a clear evidence of the approach 
to a magnetic phase transition as the temperature is de- 
creased. This critical behavior at T ~ 2.5K indeed cor- 
responds to the temperature of the cusp in the magnetic 
susceptibility and to the specific heat anomaly. It can 
therefore be associated to the occurrence of the spin-glass 
phase (see Figs. |l| and ^ in CuGa204. 
For spin-glass systems, the stretched exponential form for 
the electronic-spin contribution of the niupn-spin depo- 
larisation function Gz (t) has been shownEil to match the 
Kohlrausch-like stretched exponential for the local mo- 
ments autocorrelation function itself, which in turn arises 
from a broad distribution of electronic-spin correlation 
times. In the particular case of moderately concentrated 
systems, the exponent f3 reaches the value of | at Tf. 
On the other hand for conventional magnetic systems, 
the dynamic muon-spin depolarisation function assumes 
an exponential form (i.e. /3 = 1), reflecting the unique 
spin-relaxation frequency of the localized moments. The 
situation observed here for CuGa2 04 appears somewhat 
intermediate with an exponent /3 slightly, but definitively, 
below unity (— |). This behavior is tentatively ascribed 
to the high concentration of local moments (Cu2+ ions), 
randomly distributed in different sublattices, for which a 
somewhat narrow distribution of electronic-spin correla- 
tion times could be expected. 

In the temperature range between 3.8 and 10 K, the best 
fits with Eq. |^ provide a parameter A„s for the KT func- 
tion (i.e. essentially the width of the internal fields aris- 
ing from the nuclear moments) which is practically con- 
stant, indicating that the nuclear moments remain static 
at all temperatures. This is also confirmed by measure- 
ments performed in applied longitudinal fields (LF). If 
the nuclear moments are static within the /iSR time win- 
dow and if the applied field is sufhciently strong to quench 
the nuclear dipole field contribution (i.e. Gxrit) = 1), 
the muon-spin depolarisation should arise solely from the 
dynamical electronic-spin contribution and the depolar- 
isation function will assume the form Gz(t) = Gesit). 
This was indeed observed during LF measurements (see 
Fig. |ri|) for which a magnetic field of 0.02 T was suffi- 
cient to quench the nuclear dipolar moments. The muon 
depolarisation function is then well reproduced with the 
stretched exponential function described before, with pa- 
rameters compatible with the ones extracted from the 
zero-field data. 

For temperatures below T — 3.8 K, the muon depo- 
larization increases significantly and assumes a Gaussian 
character at short times. For this temperature range, the 
best description of the data is obtained using the function 
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ment fluctuates at the rate v. The first term 

of Eq. is only present in the temperature interval 



GS) = A para 

G GoKTit) , (10) 

where Gz, par ait) is defined above and GoKrii) ^the so- 
called dynamical Kubo-Toyabe (DKT) functionE£l, which 
refiects that the Gaussian internal field distribution due 
to the occurrence of static electronic-spins (second mo- 

between 3.8 and 2.5 K, i.e. in a region where paramag- 
netic domains appear to coexist with domains exhibit- 
ing static, albeit disordered, magnetic moments. Fig- 
ure ^ shows the temperature evolution of the amplitude 
Amagn which mirrors the volume of the magnetic do- 
mains. Therefore, it appears that in CuGa204 the tran- 
sition to a spin-glass state begins around T ~ 3.8 K to 
form local clusters of frozen electronic spins which grow 
when the temperature is lowered and finally percolate 
at the same temperature where the specific anomaly is 
observed (i.e. T ~ 2.5 K) and which can be therefore 
associated to Tf. 

With the exception of some limiting cases, the DKT func- 
tion cannot be expressed analytically and depends di- 
rectly on the parameters v and Ag^. Figure ^ shows 
the temperature dependence of the parameter Aes which 
exhibits a clear increase below ~ 3.5 K and can be associ- 
ated to the temperature dependence of the static part of 
the electronic magnetic moments. The fluctuation rate 
V was found to be constant below Tf [v ~ 3.7 MHz). It 
is worthwhile to note that the DKT function, which ap- 
pears to describe perfectly the data for T <C T/, assumes 
a single fluctuation rate v for the internal fields sensed by 
the muon spin. This has to be connected to our simple 
picture that the slightly reduced value in the paramag- 
netic phase of the exponent /3 compared to unity must 
be related to a rather narrow distribution of electronic- 
spin correlation times. Interestingly, the DKT function 
describes the data more satisfactorily than the model of 
"coexistingpStatic and dynamical fields" developed by Ue- 
mura et wj.tZl, and based on the theory of Edwards and 
AndersonQ, where each local moment is taken as a vector 
sum of a static component and a dynamical component 
randomly fiuctuating (see Fig.p^. This indicates that in 
CuGa204, by decreasing the temperature, an increasing 
part of each Cu^+ moment become quasi-static (charac- 
terized by a slow fluctuation rate v) , whereas the remain- 
ing part does not affect the muon polarization due to fast 
fluctuations that are not accessible within the /j,SR time 
window. 



VI. MONTE-CARLO SIMULATIONS 

The experimental observations presented in the pre- 
ceding sections all indicate that the Cu^+ magnetic mo- 
ments in CuGa204 undergo a phase transition to a spin- 
glass state at Tf=2.5 K. To understand the nature of 
this magnetic state, Monte-Carlo simulations were per- 
formed using a model of Heisenberg spins with compet- 



ing exchange interactions including random anisotropics. 
These arise as the result of Jahn- Teller distortions 
of the octahedcons and tetrahedrons surrounding the 
Cu^+ positionsEZl. The local distortions occur randomly 
along one of the three equivalent C4 cubic axes. Con- 
sequently, the exchange interactions between nearest- 
neighbors spins located on tetrahedral (A-sites) and oc- 
tahedral (B-sites) positions have tetragonal anisotropy. 
However, the direction of the tetragonal axis is ran- 
dom in a crystal with cubic symmetry. For the model 
calculations, we considered exchange interactions be- 
tween nearest-neighbors Cu^+ (A) — Cu^"*" (B) and second- 
nearest- neighbors Cu^+(B) — Cu^"'"(B) magnetic ions. 
We took into account the fact that in the spinel lattice the 
Cu^"*" ions are randomly distributed between the A and 
B sites with occupation probabilities of 25% and 75%, 
respectively. Consequently, the model Hamiltonian for 
this spin system is given by 



TT \ ^ \ ^ jaa na Qa pi pO 

a 2^ '-'i] '^3 

OL—x.y,z i,j 

- ^ K,,S,S,PlP^ - ^ HStiPfPf), (11) 

i,i i 

where the J^"" 's represent the exchange integrals between 
the nearest-neighbors Cu^+ ions located in the A and B 
sites; Kij is the exchange parameter between nearest- 
neighbors Cu^+ ions on the octahedral sublattice and 
H the external magnetic fleld. The components of the 
exchange interactions J,"" are distributed randomly with 
the same probability, namely 



p( jxx ryy jzz\ __ 

i/mj-; -Jo- Amjfy Jo)s{j:f jo) + 

i/36{j^; - Jo)d{j!y -Jo- Aj)S{j:'f - Jo) + 

AJ), 



l/35(Ji5-- - Jo)S{J!J' - Jo)S{J, 



yy 

ij 



Jq 



(12) 



where S{x) is the (5- function. The random numbers P/ 
and Pj determine the distribution P of the Cu^+ ions 
among the tetrahedral and octahedral sites in the spinel 
lattice, respectively, so that 

P(p;'") = iy''°S{P-'° - 1) + (1 - iy''°)6{P-'°) (13) 

with J/* = 0.25 and t/" = 0.75. The Monte-Carlo simu- 
lations were carried out using periodic boundary condi- 
tions for a lattice consisting of 24 x 24 x 24 sites and over 
30000-60000 MK steps per spins. We calculated the mag- 
netization of the spin lattice, the magnetic susceptibility, 
and the spin-spin correlation function (S(O)S(R)). We 
simulated the temperature dependence of the EA-order 
parameter q(T) for the A- and B-spins, respectively, de- 
flned as 



a/3 



(l/7V)^(Sr)2,a = x,y,z,/3 = i,0. (14) 



The exchange parameters A J, K and Jq were obtained 
by fitting the MK results to the experimental freezing 



5 



temperature Tf, the paramagnetic Neel temperature Q 
and the magnetic field dependence of the magnetiza- 
tion M(H). Fig. |l^ shows the temperature dependence 
of the magnetic susceptibility calculated by the Monte- 
Carlo method for two values of magnetic fields, H=0 Oe 
and H=10'* Oe, respectively. In agreement with the ex- 
perimental results, the calculated magnetic susceptibil- 
ity exhibits a cusp at Tf ~ 2.5K which is suppressed 
when a magnetic field is applied. As shown in Figjl^ the 
EA-order parameters for both the A- and B-sublattice 
sharply increase below T=Tf. Moreover, the spin-spin 
correlation function (S(0)S(L/2)) (L=Monte-Carlo sam- 
ple size) reveals the absence of any long-range magnetic 
ordering in the spin system. The Monte-Carlo results 
show that for the concentrations of Cu^+ spins of rele- 
vance for CuGa204, the crystal is in a superparamagnetic 
state when K=0. Introducing random anisotropy for the 
exchange interaction J results in a spin cluster block- 
ing at the freezing temperature Tf. In that respect, we 
note that the freezing temperature is proportional to A J. 
To reproduce the magnetization data in a satisfactory 
way, we found necessary to give a non-zero value to the 
antiferromagnetic exchange interaction K. From a least- 
square refinement of the field dependence of the magne- 
tization we obtained the parameter values AJ/Jo=0.1, 
K/ Jo~0.5. The exchange parameter Jo, as determined 
from the freezing temperature T/ , from the paramagnetic 
susceptibility in the temperature range 90K < T < 160K 
and from the magnetization curve, amounts to -12K, - 
12. 5K and -13 K, respectively. Therefore, the mean val- 
ues of the model parameters are J=-12.5K, AJ=-1.3K, 
and K=-6.2K. 



VII. CONCLUSION 

We have presented magnetization, magnetic suscepti- 
bility, specific heat and /iSR measurements in CuGa204. 
The data are consistent with a spin-glass transition of the 
copper sublattice below Tf in this material. In particu- 
lar, we observe a cusp in the temperature dependence of 
the magnetic susceptibility at Tf ~ 2.5K which is sup- 
pressed when a magnetic field is applied. A pronounced 
hysteresis is observed in the temperature dependence of 
the magnetic susceptibility for zero-field and field-cooled 
samples. The muon-spin relaxation measurements have 
shown that above the freezing temperature, the asym- 
metry function is described by the stretched exponential 
typical of disordered systems. However, the value of the 
exponent f3 points to a narrow distribution of correla- 
tion times of the local moments. The temperature de- 
pendence of the magnetic volume fraction indicates that 
in CuGa204 the transition to a spin-glass state begins 
around T~3.8K to form locally clusters of frozen spins 
which grow when the temperature is lowered and finally 
percolate around Tf ~ 2.5 K. By means of Monte-Carlo 
simulations we were able to reproduce the main features 



of the magnetic susceptibility and of the magnetization 
curve measured in CuGa2 04. The results of Monte-Carlo 
simulations show that the Jahn- Teller effect plays an es- 
sential role in forming the magnetic ground-state as it 
introduces random anisotropy in the exchange interac- 
tions between the copper ions, which in turn is responsi- 
ble for the formation of a spin-glass state in CuGa204. 
Using a realistic spin model which takes into account the 
effective distribution of the Cu^"'' ions in the spinel struc- 
ture, reliable exchange parameters could be obtained for 
CuGa204. 
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FIG. 1. (a) Real and (b) imaginary parts of the magnetic 
susceptibility in a single crystal of CuGa2 04 as measured with 
a SQUID magnetometer at a frequency of 19 Hz and in an 
applied magnetic field of 4.5 Oe. 



FIG. 10. Temperature dependence of the depolarisation 
rate A above 3.8 K. 

FIG. 11. ^SR spectra measured in CuGa204 at T = 5 K 
and showing the field dependence of the asymmetry function. 
Whereas in zero-field the depolarization function is best de- 
scribed according to Eq. ^ the data with applied field are 
only fitted by the stretched exponential term refiecting the de- 
polarisation arising from the fiuctuating electronic spins (see 
text). 

FIG. 12. Temperature dependence of the parameter Amagn 
corresponding to the magnetic volume fraction. 

FIG. 13. Temperature dependence of the A-parameter of 
the DKT function. This parameter mirrors the width of the 
quasi-static field distribution below Tf and therefore the value 
of the quasi-static Cu^^ moment. 



FIG. 2. Experimental value of the magnetic moment in 
CuGa2 04 for a sample cooled in an applied magnetic field 
of H=100 Oe (FC curve) and in zero-field (ZFC curve). 



FIG. iM. Fit to a ^SR spectrum using the Uemura's 
functioncJ (see also text). For sake of clarity the error bars 
are omitted. The shortcomings of the fit are evident. 



FIG. 3. Temperature dependence of the magnetic suscep- 
tibility for different magnetic fields. 

FIG. 4. Curve 1 shows the magnetic field dependence of the 
magnetic moment in CuGa204 at the temperature of T=1.8 
K. The magnetic field is applied along the [001] crystal axis. 
In curve 2, we show the result of the Monte-Carlo simulation 
for the same field dependence. 

FIG. 5. Magnetic susceptibility measured in CuGa204 for 
different frequencies. 



FIG. 15. Temperature dependence of the normalized mag- 
netic susceptibility for values of magnetic fields H=0 Oe and 
11=10* Oe (curve 1 and 2, respectively) with exchange con- 
stants J=-12 K, K=-6 K and AJ=-1.2 K. 

FIG. 16. Temperature de- 

pendence of the Edwards- Anderson parameter q(T) for the 
A- and B-sublattice as simulated with Monte-Carlo. See text 
for details. 



FIG. 6. Specific heat of CuGa204. 



FIG. 7. Temperature dependence of M/H in CuGa204 at 
T=1.8 K. 



FIG. 8. Experimental and theoretical dependences of the 
freezing temperature Tf as a function of the reduced magnetic 
field h, as explained in the text. 



FIG. 9. Experimental zero-field /iSR signal measured in 
CuGa204 at T = 10 K, T = 4.5 K and T = 650 mK. The 
lines represent fits as explained in the text. 
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